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ABSTRACT 

Recent spectroscopic observations have revealed the ubiquitous presence of 
blueward asymmetries of emission hues formed in the solar corona and transi- 
tion region. These asymmetries are most prominent in loop footpoint regions, 
where a clear correlation of the asymmetry with the Doppler shift and line width 
determined from the single Gaussian fit is found. Such asymmetries suggest at 
least two emission components: a primary component accounting for the back- 
ground emission and a secondary component associated with high-speed upfiows. 
The latter has been proposed to play a vital role in the coronal heating process 
and there is no agreement on its properties. Here we slightly modify the ini- 
tially developed technique of Red-Blue (RB) asymmetry analysis and apply it 
to both artificial spectra and spectra observed by the EUV Imaging Spectrom- 
eter onboard Hinode, and demonstrate that the secondary component usually 
contributes a few percent of the total emission, has a velocity ranging from 50 
to 150 km s~^ and a Gaussian width comparable to that of the primary one in 
loop footpoint regions. The results of the RB asymmetry analysis are then used 
to guide a double Gaussian fit and we find that the obtained properties of the 
secondary component are generally consistent with those obtained from the RB 
asymmetry analysis. Through a comparison of the location, relative intensity, 
and velocity distribution of the blueward secondary component with the prop- 
erties of the upward propagating disturbances revealed in simultaneous images 
from the Atmospheric Imaging Assembly onboard the Solar Dynamics Observa- 
tory, we find a clear association of the secondary component with the propagating 
disturbances. 
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Introduction 



EUV spectroscc 
aries of some active 


)pic observations often reveal pate 


les of b 


dueshifted emission at bound- 


regions (ARs, '. 


Vlarsch et al.l 2004 


. 2008 


Harra et al. 2008: Del Zanna 


2008: Dosche 


c et al. 


2008: Triuathi et al. 


2009: Murrav et al. 


2010l: Brooks & WarrenI 2011: 


Warren et al. 


2011: 


Del Zanna et al. 


2011 


). These authors used a single Gaussian fit to ap- 



proximate coronal emission line profiles and derived a blue shift of the order of 20 km s 



Large line widths have also been found in these blueshifted regions (e.g., iDoschek et al. 



20081 ) . There are suggestions that these b l ue shifts are indicators of the nascent slow s olar 
wind dSakao et allboO?! : iHarra et al.lbood : IPoschek et aPbood : brooks fc WarrenI boilh . 



Recent investigations have revealed the ubiquitous presence of blue ward asymmetries 



of emission lines formed in the solar corona and transition region (e.g., iDe Pontieu et al. 



20091 : iMcIntosh fc De Pontieul l2009bl ) . The enhancement in the blue win gs of line profiles 



are most prominent in loop footpoint regions, such as AR boundaries (IHara et al 



De Pontieu et al 



2009 



De Pontieu fc Mcintosh 



2010 



Mcintosh fc De Pontieu 



2009a 



20081: 



Peter 



2011al :l Ugarte-Urra fc Warre J [2OI1I ) . The discovery of 



2OIOI : iBrvans et al. 



I2OI 



Tian et al. 



these asymmetric line profiles suggests the presence of a highly blue-shifted (much larger 
than 20 km s~^) emission component besides the primary emission component, and thus 
provides a significant challenge to previous results and interpretations which are based on a 
single Gaussian fit. 



De Pontieu et al.l ( 120091 ) and iMcIntosh fc De Pontieul ( j2009bl) suggested that this sec 



ondary emission component is a ssociated with type-II spicules ( 



rapid blue-shifted events (RBEs, iRouppe van der Voort et al. 



De 



2009 



Pontieu et al 



2007h or 



De Pontieu et al 



2OIID 



plasma (e.g.. iDe Pontieu et al.l l2009l : IMcIntosh fc De Pontic 



De Pontieu et al. 



a rep 
bieul I 



Hansteen et al. 


2010 


^sis ( 


De Pontieu et al. 



2OO9I : iMartmez-Sykora et al.l l2011al ) is based on a comparison of the two wings at same 
velocity ranges, and detailed analysis indicates that both the speed (50-150 km s~^) and 
the relative intensity (a few percent of the total emission) of the secondary component 
are roughly consistent across a temperature range from 100,000 to several million degrees 



( IDe Pontieu et al. 



2009 



Mcintosh fc De Pontieul l2009br> . Through j oint imaging and spec- 



troscopic observations of the corona. IMcIntosh fc De Pontieul (j2009al ) and lTian et al.l ( l2011al ) 
have suggested that the secondary emission component is caused by high-speed repeti- 
tive upfiows in the form of upward propagating disturbances in EUV and X-ray imag- 
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ing observations. S uch disturbances were previously interpreted as slow mode maRnetoa- 



coustic waves (e.g.. iDe Moortel et al. 



^ere p 
200oi 



2002 



Robbrecht et al 



2001 



McEwan fc De MoortefcoodlMarsh et al.lbood i lWang et al.ll2009Jbl:lDe Moortel 



King et al 



2003 



2011 



20101 : 



Marsh et al 



Guo et al 



2011 ) or slow-speed solar wind outflows (e.g., Sakao et al. 200?! : Ihc et al. 



2009; 



Stenborg et al. 



2010[>. There is also a suggestion that they are warps in two-dimensional 



sheet-like structures ( 1 Judge et al.ll201ll ). 



Besides the RB asymmetry analysis, the double Gaussian flt technique has also been 
used to resolve the second emission component from li ne profiles (IDe Pontieu fc Mcintosh 



20101: 



Pete: 



2010l:lBrvans et al.lboiol : Ixian et al.l[2011al ). be Pontieu fc Mclntoshl koid ) and 



Tian et al.l ( l2011al ) applied a RB guided double Gaussian fit to spectra with a high signal- 
to-noise ratio and obvious blueward asymmetry. They used the velocity derived from the 
RB asymmetry analysis as an initial guess of the velocity of the secondary component. The 
algorithm undertakes a global minimization of the difference between the observed spectrum 
and the fit by allowing both primary and secondary component centroids to move by one 
spectral pixel (^30 km s~^) to the blue or red of the initial positions. The s p eed o f the 



secondary com ponent was found to be ^60 km s ^ in iDe Pontieu fc Mclntoshl ( 120101 ) and 



100 km s ^ in lTian et al.l ( l2011al ). They both found that the widths of the two components 



are similar. The widths of the tw o Gaussian co r apone nts were forced to be the same in the 
double Gaussian fit algorithm of i Bryans et al. ( 2010 ). who mentioned that the velocity of 
the secondary component is often as high as 200 km s~^ and that the primary component 
is also blue-shifted by ^10 km s~^. However, they claimed that the double Gaussian fit is 
a good approximation for only the Fe xil 195. 12A and Fe xill 202. 02A lines and that other 
lines are fine with a single Gaussian fit. By app l ying a completely free double Gaussian fit to 
spectral profiles of the Fe XV 284A line, IPeterl ( 120101 ) found that the secondary component 
usually contributes 10% to 20% to the total emission, is usually blue-shifted by ^40 km s~^ 
and twice as broad as the primary component, especially in loop footpoint regions. However, 
different settings of initial values and allowable ranges of the seven free parameters can lead 
to very different fitting results, and thus reasonable initial values and constraints to some 
free parameters are highly desirable when performing the double Gaussian fit. 

Here we slightly modify the initially developed technique of RB asymmetry analysis 
and apply it to both artificial and observed spectra by the EUV Imaging Spectrometer (EIS, 
Culhane et al.ll2007l ) onboard Hinode^ and demonstrate that the speed of the secondary com- 
ponent usually does not reach 200 km s~^ and that the widths of the two components are 
comparable. We then use parameters determined from the RB asymmetry analysis to guide 
the double Gaussian fit to each spectrum that is observed to have an obvious blueward asym- 
metry. Both the RB asymmetry analysis and the double Gaussian fit yield consistent results. 
We also use imaging observations simultaneously performed by the Atmospheric Imaging 
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Assembly (AIA, iBoerner et al.ll2010l ) onboard the Solar Dynamics Observatory (SDO) to 
demonstrate that the propagating disturbances at AR edges are plasma upflows which are 
indeed responsible for the blueward asymmetries of line profiles. 



2. Artificial line profiles 



The technique of RB a symmetry analysis was first introduced bv iDe Pontieu et al 



(12009). In our previous work (IDe Pontieu et al 



2009: 



De Pontieu fc Mcintosh 



3 



201 



De Pontieu et al 



2011 



Mcintosh fc De Pontieu 2009a 



Mcintosh et al 



2010 



Tian et al.ll2011al : lMartmez-Sykora et al 



2011al ). we first interpolated the line profile to a spectral resolution ten times greater than 
the original one, then subtracted the blue wing emission integrated over a narrow spectral 
range from that at the same position and over the same range in the red wing. The range 
of integration was then sequentially stepped outward from the line centroid to build an RB 
asymmetry profile (simply RB profile). The RB a s ymmetry for an offset velocity Uc can be 



expressed as the following (IDe Pontieu et al.l l2009l : iMartmez-Sykora et al.l l2011al ) : 



'uc—Su/2 



-Uc+Su/2 
Uc—Su/2 



I{u)du^ 



(1) 



where 5u^ and I{u) represent the velocity from the line centroid, the velocity range 
over which the RB asymmetry is determined, and the spectral intensity, respectively. The 
value of 8u is usually set as 20 km s~-^. In our previous work, we used the single Gaussian 
fit to determine the line centroid. The resulting RB profile can be normalized to the peak 
intensity derived from the single Gaussian fit. 

Here we slightly modify this technique by using the spectral position corresponding 
to the peak intensity as the line centroid. The resulting RB profile is then normalized to 
the peak intensity of the observed line profile. In the following the RB profiles all refer to 
normalized RB profiles. The originally defined and modified techniques are designated as 
RB5' and RBp, respectively. In the following we apply these two techniques to artificial 
spectra composed of two Gaussian components to test the ability of the RB asymmetry 
analysis to resolve the secondary component. 



Similar to iMartmez-Sykora et al.l (l2011al ) , each artificial spectrum consists of a fixed 
primary component and a much weaker secondary component. However, to mimic the ob- 
served EIS spectra, th e Gauss ian parameters we used here are very different from those in 
Martinez- Sykor a et al.l (j2011al ). The spectral pixel size is set as 24 km s~^, similar to that 
of the EIS spectra at around 274A. The relative peak intensity of the secondary component 
with respect to the primary one can be set as any value below 30%, while in Figured] we only 
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Fig. 1. — Parameters derived from the RB asymmetry profiles of artificial emission line pro- 
files. First row: the primary, secondary, and total emission profiles are shown as the dashed, 
dash-dot, and solid lines, respectively. The peak intensity of the secondary component is set 
as 20% of that of the primary one. Second row: RB^- (solid) and RBp (dashed) asymme- 
try profiles. Different colors represent different velocity offsets of the secondary component. 
Third row: the relative intensity of the peak of the asymmetry profile. Fourth row: the 
velocity and 1/e width of the asymmetry profile. Solid and dashed lines represent results of 
RB^- and RBp, respectively. The Gaussian width of the secondary component is 28, 56, and 
112 km s~^ for the left, middle, and right columns, respectively. 
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show the case of 20%. The Gaussian width (1/e width) of the primary component is set to 
be 56 km s~^, a value comparable to the coronal line width (including instrumental, thermal, 
and non-thermal width) observed by EIS. For the secondary component, three values of the 
width are chosen (28, 56, 112 km s~^). The spectral position of the secondary component is 
shifted by -20, -40, -60, -180 km s~^ with respect to the major component and they are 
shown as different colors in Figure [H 

After applying the RB^- and RBp techniques to each artificial spectrum, we obtain RB 
asymmetry profiles of each spectrum. The same colors are used in the RB asymmetry plots 
of Figure [H In the third and fourth rows of Figure [1] we plot the peak intensity, velocity and 
1 /e width derived from the RB analysis as a function of the velocity oflPset of the secondary 
component. 



As noticed by iMartmez-Sykora et al.l (l2011al ) , when the width of the secondary compo- 



nent is considerably smaller (28 km s~^) than that of the primary one and the secondary 
component is centered at velocities lower than the Gaussian width of the primary component, 
the RB^- asymmetry profile is blueward at low velocities and redward at higher velocities. 
The asymmetry profile reverses when the offset velocity of the secondary component is larger 
than the width of the primary component. Here we find that the RBp asymmetry profile is 
redward at low velocities and blueward at higher velocities when the secondary component 
is too close to the primary one. As the offset velocity increases to a value larger than the 
width of the primary component, the RBp asymmetry profile is blueward at all velocities. 
We find that when the velocity offset is larger than the width of the primary component, 
the velocity, width, and relative intensity derived from the RBp technique are very close to 
the true velocity and width of the secondary component. The RB^- technique does slightly 
worse but still can give relatively accurate values of the parameters. 

When the widths of two components are the same (56 km s~^), no matter where the 
secondary component is centered, the RB^- asymmetry profile is redward at low velocities and 
blueward at higher velocities, and the RBp asymmetry profile is blueward at all velocities. 
By comparing the velocity, width, and relative intensity derived from the RB techniques with 
the true parameters of the secondary component, we can still conclude that the RB tech- 
niques can provide relatively accurate estimate of the properties of the secondary component 
provided the offset velocity is larger than the width of the primary component. We notice 
that compared to the RB^- technique, the RBp technique has a better ability to reproduce 
the parameters of the secondary component, especially when the offset velocity is smaller 
than twice the width of the primary component. The improvement is around 10 km s~^ in 
the velocity and width, and 4% in the relative intensity. 

When the secondary component is considerably broader (112 km s~^) than that of the 
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primary one, the RB profiles show morphologies similar to the case of 56 km s"-*^. However, 
both the RB^- and RBp techniques can not reproduce an offset velocity smaller than 100 
km s~^. The velocity derived from RB analysis is always larger than 100 km s~^ if the 
secondary component is twice as broad as the primary one. There is also a large discrepancy 
between the width derived from RB analysis and the true width of the secondary component 
if the offset velocity is smaller than twice the width of the primary component (112 km s~^). 

We notice that in all the three cases, the RB techniques can give an extremely accurate 
estimate of the velocity when the offset velocity of the secondary component is larger than 
112 km s-^ 

We have varied the intensity ratio of the two components in the range of 5-30% and 
obtain basically the same results, except for the different values of the calculated relative 
intensities. One thing we have noticed from Figure [1] is that the RB technique tends to 
underestimate the value of the relative intensity of the secondary component at smaller 
offset velocities (e.g., 56-112 km s~^). This is especially the case for the RB^- technique 
which we used in our previous work. From the third row of Figure [1] we can see that in the 
velocity range of ^56-112 km s~^ the relative intensity of the secondary component recovered 
from the RBp technique is much closer to the real value (-0.20 in Figured]), as compared to 
the RB^- technique. Th e underestimation of relative intensity might be improved by using 



the Chi-square method (IMartmez-Sykora et al.ll2011al ). 



The fact that the calculated RB velocities and relative intensities obviously deviate from 
the real values when the offset velocities of the secondary component are smaller than the 
width of the primary component has an important implication for the development of future 
EUV spectrographs. The EIS instrument has a large instrumental width of ^35 km s~^ (e.g.. 



Doschek et al.l 120071 ). which contributes almost half of the total line width. The Interface 
Region Imaging Spectrograph (IRIS), which is scheduled to be launched at the end of 2012, 
will have a much smaller instrumental width (a Gaussian width of ^5 km s~^ at 1400 A) and 
our RB technique should be able to accurately resolve the secondary component at smaller 
offset velocities. The very high spectral resolution (^3 km s~^ ) of the IRIS spectrograph 
will greatly reduce the error caused by the interpolation of line profiles. 

In the following we pick out two ARs and investigate properties of the secondary com- 
ponent from the spectra obtained by EIS. Since the modified RB technique, RBp, has a 
better ability to accurately resolve the blueshifted secondary component as compared to 
the originally defined RB^- technique, here we mainly present results by appling the RBp 
technique instead of the RB^- technique to the real data. A comparison between the RBp 
and RB^- results for the 2010 September 16 observation (see the details of this observation 
in Section 4) is presented in Appendix A. 
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3. EIS observation on 2007 January 18 



The I" X 512^' slit of EIS was used to scan AR NOAA 10938 from 18:12 to 20:27 on 
2007 January 18, with an exposu r e tim e of 30 s and a step size of 1" . This data was 
previously analyzed by iHara et al.l (120081) andlPeter for the purpose of investigating 

properties of asymmetric line profiles. Hara et al. ( 2008 1) mainly used the Fe xiv 2 7 4. 20 A 
and Fe XV 284. 16 A lines which are formed around a temperature of 2 MK, while IPeterl ( 120101 ) 
only focused on the Fe X V 284. 1 6A lin e. Here we choose both lines for our study but focus on 
the Fe xiv 274. 20A line. IPeterl (I2OI0I ) mentioned that the Fe XV 284. 16A line is blended with 
Al IX 284. 03A and that the latter usually contributes no more than 5% to the deviation from 
a single Gaussian profile. However, as we will discuss later, in some locations the blend can 
greatly complicate the asymmetries of line profiles. The Fe xiv 274. 20A line is also a strong 
line and it is blended with Si Vll 274. 18A. The blend is much weaker than Fe xiv 274. 20A 



and can safely be ignored in active region conditions (l Young et al.l 120071 ). Moreover, this 
blend is very close (less than 1 spectral pixel) to the line center of Fe xiv 274. 20A and thus 
it should not have an important infiuence on the results of our RB asymmetry analysis and 
Gaussian fit. Further confidence is given by the highly similar behavior of the line moments 
and profile asymmetries between Fe xiv 274. 20A and the weaker Fe xiv 264. 78A line in most 
ARs we analyzed (see Appendix B), although the Fe xiv 264. 78A line is usually noisier than 
Fe XIV 274.20A. 

The SSW routine eis_prep.pro was applied to correct and calibrate the EIS data. This 
includes CCD pedestal and dark current subtraction, cosmic ray removing, warm and hot pix- 
els identification, absolute calibration, error estimation, and so on. The eflFects of slit tilt and 
orbital variation (thermal drift) were estimated by using the SSW routine eis.wave.corr.pro 
and removed from the data. After that, a running average over 5 pixels along the slit and 3 
pixels across the slit was applied to the spectra to improve the signal to noise ratio. 

As a common practice, a single Gaussian fit was applied to each spectrum and Fig- 
ure [2](A)-(C) show the spatial distributions of the Gaussian parameters for Fe xiv 274. 20A. 
Figure [2](D) s hows the map of reduced xl for the fit. The reduced Xr defined as the 
following (e.g.. iBevington fc R.obinsonlll992l : IPeterl 120011 ): 



Xr-l^N^ — ^ — ' y^) 

where d^, m^, and denote the observed spectral radiance, fitted spectral radiance, and 
measurement error calculated by using the SSW routine eis_prep.pro (mainly Poisson error), 
respectively. Here i represents the spectral position and the summation is performed over 
all the spectral positions. The degree of freedom is given by / and equals to 4 or 7 for a 
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(A) Peak Intensity (arbitrary unit) (^) Doppler Shift (km/s) (^) Gaussian Width (km/s) (^) of Single Gaussian Fit 
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Fig. 2. — Spatial distributions of the parameters derived from the single Gaussian fit (A-C) 
and RB asymmetry analysis (RBp, E-G) for Fe XIV 274.20A. The map of Xr is shown in 
(D). The map of the average relative intensity in the velocity interval of 60-110 km s~^, as 
obtained from the RB profiles, is shown in (H). The two squares in (A) mark the locations 
where profiles are averaged and presented in Figure [H 



-10- 



single or double Gaussian fit. 

We assume zero shift of the profile averaged over the entire observation region. We can 
see that the loop footpoint regions, or boundaries of the AR, are characterized by a blueshift 



Hinode era (e.s.. Marsch et al. 


2008: Harra el 


: al. 20081: Del Zanna 2008: Doschek et al. 2008: 


Mcintosh & De Pontieul 


2009a 


:lMurrav et al. 


2010l: Peter 201ol:lTian et al. 201 lal: Brooks & Warren 


2011: Warren et all 


2011 


: Del Zanna et al. 201lh. However, the goodness of the fit. the Xr 



is also clearly enhanced mainly in loop f ootpoint regions, indicating an obvious deviation 



fr om a s i ngle G aussian profile there (e.g., iPeter 



m 



2001 



2010|). Note that the definition of x: 



Peterl ( 120101 ) is not correct and that it is probably a typographic error. 



3.1. RB asymmetry analysis 

We applied the modified RB technique RBp to each spectrum and calculated the average 
relative intensity in the velocity interval of 60-110 km s~^ from the RB asymmetry profile. 
Figure [2](H) shows the spatial distribution of this average RB asymmetry. Here a nega- 
tive/positive value indicates an enhancement of the blue/red wing in this velocity interval. 
We can see that the blueward asymmetry is most prominent at loop footpoint regions, gen- 
erally coincident with the blueshift, width enhancement, and Xr enhancement. For further 
analysis, we only selected those locations where this average RB asymmetry is smaller than 
-0.05 (obvious blueward asymmetry) and the signal to noise ratio of the profile (defined as 
the ratio of the peak and background intensities) is larger than 5. We have to mention that 
such a cutoflF excludes the analysis of spectra in the AR core and edges where the magnetic 
field lines largely incline with respect to the line of sight, since an inclination of the magnetic 
field line usually leads to a smaller offset velocity of the secondary component in the line of 
sight direction, which would reduce the value of the relative intensity derived from the RB 
technique (see Figured]). Figure [2](E)-(G) show the peak relative intensity, velocity, and 1/e 
width derived from RBp asymmetry profiles at these locations. The distributions of these 
three parameters are presented as the red histograms in Figure [6]( A)- (C). 

In Figure [3] we show the scatter plots of the relationship between Doppler shift /Gaussian 
width derived from the single Gaussian fit and the peak relative intensity derived from RBp 
asymmetry profiles for Fe XIV 274. 20A. Here only data points with blue shift at those se- 
lected locations are shown. A clear correlation is found in each of the two panels in Figure [3l 
Such a correlation strongly suggests that the clear blue shift and enhanced line width in 
the loop footpoint regions are related to, or caused by the blueward asymmetries. We can 
simply imagine a faint high-speed upfiow superimposed on a strong and almost stationary 



- 11 - 



(or slightly shifted) background in the line of sight. Such a scenario would naturally lead 
to a blueward asymmetric line profile. And a single Gaussian fit to the total emission line 
profile would result in a blueshift and enhanced line width, as com pared to the line pro- 
file of the background emission-a similar conclusion was reached by [PeterJ ( 120101 ). So it is 



clear th at the blue shift at AR boundaries, which was reported in many previous investi- 



Rations ( 


Marsch et al. 


2004, 


2008; 


Harra et al. 


2008; 


Del Zanna 


2008; 


Doschek et al. 


2008; 


Trioathi et al. 


2009: A 


/lurrav et al. 


201 


0; 


Peterl 


2010; 


Brooks & Warren 


2011; 


Warren et al. 


2011 


: lOel Zanna et al 


2011; 


He et al. 


2010 


), is actually a composite effect of at least two 



emissi o n components and can not refl e ct the real physical process (IMcIntosh fc De Pontieu 



2009al : 



De Pontieu fc Mcintosh 



20101 : 



Peter 



real p . 
201ol : 



Brvans et all l2010l : IXian et all l2011ah . 



Some authors claim that the blue shifts derived from the single Gaussian flt at AR bound- 



aries are indicators of the nascent slow solar wind (ISakao et al.l 120071 : iHarra et al.l 12008 



Doschek et al.ll2008l ). If this is the case, the wind speed must have been largely underesti- 
mated if the source of the mass supply to the solar wind is the secondary component, or 
overestimated if the primary component dominates the mass supply to the solar wind. Our 
results also suggest that the enhanced line width can be the result of the superposition of 
different emission components with di fferent velocities and i s not necessarily related to the 
increase in the Alfven wave amplitude ( IDoUa fc Zhukovll201ll ). We think that both the inho- 
mogeneities of flow velocities and Alfven waves may contribute to the observed non-thermal 
width of an emission line, and that caution m ust be taken when interpreting the non- thermal 
width as purely cased by Alfven waves (e.g., iBanerjee et aPllQQsl . boool : IPeteilboioh . 



From Figure [6]( A) we can see that the distribution of the intensity ratio (relative intensity 
of the secondary component) peaks around 0.07 but can reach as high as 0.3. This is generally 
consistent with previous result that the blueward excess emission is often a few percent 
of the total emission, which is based on the RBc technique (e.g., De Pontieu et al.l l2009l : 
Mcintosh fc De Pontieul l2009bl ) . While IPeterl ( 120101 ) found a relative intensity of some 10% 
to 20% based on his free double Gaussian flt. 

As mentioned above, when the width of the secondary component is not considerably 
larger than that of the primary one, the RBp technique can provide an accurate estimate of 
the Gaussian parameters of the secondary component provided the offset velocity is larger 
than the width (~56 km s~^) of the primary component. As seen from Figure[6](B), the RBp 
velocity seldom reaches 100 km s~^ and its distribution peaks at ^75 km s~^. Such a result 
would exclude the possibility of a very broad secondary component, as claimed by IPeterl 
( 120101 ) based on a free double Gaussian flt. This is because the velocity derived from RB 
analysis is always larger than 100 km s~^ when the secondary component is twice as broad 
as the primary one (see Figured]). In Figure [6](C) we can see that the RBp width is usually 
smaller than 60 km s~^ and larger than 35 km s~^, and that its distribution peaks around 
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50 km s~^, which is also not consistent with what we expect if the secondary component is 
twice as broad as the primary one. From the last panel of Figured] we can see that the RBp 
width is always larger than 60 km s~^ in such a case. From Figure [1] we can also see that 
when the width of the secondary component is considerably smaller (28 km s~^) than that 
of the primary one, the RBp width is always smaller than 35 km s~^, which is not consistent 
with the observed histogram in Figure [6](C). These inconsistencies strongly suggest that the 
secondary component can not be considerably broader or narrower than the primary one. 
The widths of the two components should be comparable. This conclusion implies that the 
as sumption of the sam e width of the two components in the double Gaussian fit algorithm 
of iBrvans et al. is reasonable. 



The velocity distribution presented in Figure [6](B) is not consistent with the scenario o f 



a very high-speed (^200 km s~^) secondary component, as suggested by lBryans et al.l (120101 ) 
after analyzing another data set. As seen from Figured! the velocity derived from the RBp 
technique is extremely close to (usually less than 5 km s~^) the true velocity of the secondary 
component when the latter is larger than 112 km s~^. If very large velocities exist, the RBp 
technique should yield such large numbers of velocities. However, the observed velocity 
derived by using the RBp technique seldom reaches higher than 150 km s~^, indicating 
that the velocity of the secondary component in our observation can not be as high as 200 
km s~^. We have to mention that the magnitude of velocity also depends on the viewing 
angle. However, in loop foot point regions usually a large portion of the magnetic field lines 
are almost vertical so that the line of sight effect should not be very significant at disk ce nter. 



We also applied the RB asymmetry analysis to the data analyzed by lBryans et al.l ( 120101 ) and 
found larger velocities (^95 km s~^) of the secondary component in that AR (see Appendix 
C). But still the velocity seldom reaches 200 km s~^. 



3.2. Double Gaussian fit 



Following the initial idea of using; the RB.^ ana lysis result to guide double Gaussian 
fit ( jPe Pontieu fc Mclntoshl l2010l : iTian et al.l l2011al ). we designed a similar but modified 
algorithm. The peak intensity, velocity, and width derived from a single Gaussian fit were 
used as the initial guess of those parameters of the primary Gaussian component. The single 
Gaussian background was also used as the initial value of the double Gaussian background. 
Meanwhile, we used the peak intensity and velocity obtained through the RBp analysis as 
the initial values of the same parameters of the secondary Gaussian component. The width 
of the secondary component was initially set as the same as the primary component. During 
the iterations, we allowed all the three intensities (peak intensities of the two components 
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and the background intensity) to vary within the range of 75%- 125% of the corresponding 
initial values. The velocity was allowed to move to the blue or red of the initial position by 
one spectral pixel (^24 km s~^ at 274A) for the primary component and two pixels for the 
secondary component. We also allowed the width to increase or decrease by one spectral 
pixel for the primary component and two pixels for the secondary component. The algorithm 
undertook a global minimization of the difference between the observed and fitted spectrum. 

We present several examples of the observed and fitted profiles of both the Fe xiv 274. 20A 
and Fe XV 284. 16A lines in Figure [H The profiles in Figure [11(A) & and (B) are the profiles 
averaged respectively over region 1 and 2 marked in Figure [21 By comparing the observed 
profiles with the different fitting profiles, we can clearly see the better performance of the 
double Gaussian fit and the deviations of the observed profiles from the single Gaussian fits. 

After the double Gaussian fit, we took the velocity of the primary component as the line 
centroid and calculated the RB asymmetry profile (here RB^). The lower panels of Figure [H 
show the three RB asymmetry profiles (RB^-, RBp, RB^)). We can see that RBp and RBjj 
show similar behaviors and values of the three parameters are also very close. 



I t is known that the Fe xv 284. 16A line is blended with Al ix 284. 03A (e.g. jYoung et al. 



20071 : [Brown et al.ll2008l ). The latter is sitting at ^150 km s~^ away from the spectral posi- 
tion of Fe XV 284. IGA so that this blend is potentially contaminating the blue wing of the 
line and may complicate the RB analysis and Gaussian fit. For example, the line profiles in 
both Figure [11(A) and (B) were obtained around loop footpoint regions (marked in Figure [2]). 
The magnitude of asymmetry is comparable in both the Fe XIV 274. 20A and Fe XV 284. 16 A 
line profiles shown in Figure [11(A) , indicating that the blend Al ix 284. OSA can be ignored. 
However, in Figure [11(B) we see a ^2% blueward asymmetry in Fe xiv 274. 20A but a ^16% 
blueward asymmetry in Fe XV 284. IGA. Since the two lines have similar formation tempera- 
tures, the completely different magnitude of asymmetries in the two lines strongly suggests 
that the bump at the blue wing of Fe XV 284. leA is dominated by the blend Al ix 284.03A. 
There f ore, w e focus on the Fe xiv 274. 20A line instead of Fe xv 284. IGA which was used by 



Peteil torn . 



We then applied the RBp guided Gaussian fit and RB^ analysis to all spectra of 
Fe XIV 274. 20A with obvious blueward asymmetries and large signal to noise ratio. We 
present the spatial distributions of the three Gaussian parameters of both components in 
Figure [5l In addition, the intensity ratio and velocity difference of the two components are 
also shown in Figure [51 We found that in most l ocations the p r imary component is slightly 



blue shifted (^10 km s ^), a result also found b vlBrvans et al.l (l2010l ). who interpreted it as 
the nascent slow solar wind from AR boundaries ( Wang et al. 20091 ) . While for the secondary 



component, the velocity is largely blue shifted (^80 km s ^) at almost all locations. The 
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maps of parameters derived from KBjj asymmetry profiles are very similar to those derived 
from RBp and thus are not shown here. 

Histograms of the intensity ratio (relative intensity), velocity diflFerence (velocity), and 
Gaussian width of the secondary component as derived from double Gaussian fit and RB 
asymmetry analysis (RBp,RBi)) are presented in Figure [61 We can see that they are very 
similar, especially RBp and RBp). The intensity ratio peaks around 0.07 and can also reach 
as high as 0.3. For RBp and RBp), the velocity and width both have a narrow distribution 
which peaks at ^75 km s~^ and ~50 km s~^, respectively. While for the double Gaussian fit, 
most of the width values are still much smaller than ^100 km s~^ and the width distribution 
peaks slightly higher. The velocity difference is usually in the range of 50-100 km s~^, but the 
double Gaussian fit algorithm also produces some relatively small values of velocity diflFerence 
(10-50 km s~^). However, these relatively small blue shifts are found at only a few locations 
and they are not prevalent. 

The Xr value of the double Gaussian fit is usually smaller than unity-similar result was 



also found by lPeterl ( 2010 l). who attributed it to the overestimation of the measurement error. 



From Figure [6](D) we can see that the Xr ratio between the double and single Gaussian fit 
is usually less than 0.5 in the regions characteristic of obvious blueward asymmetries. Such 
a result indicates that the double Gaussian fit does at least two times better than the single 
Gaussian fit for these asymmetric line profiles. It is likely that even more components, e.g., 
a component from cooling downfiows, might be present in the observed emission. However, a 
reliable decomposition of these additional components can not be made due to the relatively 
large instrumental width and modest spectral resolution of EIS. With the upcoming IRIS 
data, we may be able to resolve more coronal emission components. 



4. EIS and AIA observations on 2010 September 16 

The other observation was performed by EIS on 2010 September 16, with simultaneous 
observations by AIA. The EIS instrument scanned NOAA AR 11106 from 10:38 to 11:57 
UT. The 2^' x 512^' slit was used for the observation, with a 60 s exposure. After standard 
correction and calibration of the EIS data, a running average over 5 pixels along the slit 
and 3 pixels across the slit was applied to the spectra to improve the signal to noise ratio. 
Four relatively strong and clean emission lines in the spectral window were selected for our 
study: Si vii 275. 35A, Si x 258. 37A, Fe xiii 202. 04A, and Fe xi v 274.20A. Their formation 



temperatures are log(T/K)=5.8, 6.1, 6.2, and 6.3, respectively ( lYoung et al.l 120071 ). But we 
will mainly concentrate on the hottest strong line, Fe XIV 274. 20A. 
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A single Gaussian fit was first applied to each spectrum to derive the line intensity, 
centroid and Gaussian width. By assuming the average Doppler shift of each line is zero over 
the entire region, we calculated the Doppler shift from the line centroid for each spectral 
profile. Maps of the Gaussian parameters and Xr 274. 20A are presented in 

Figure [71(A)- (D). Again, we clearly see that the loop footpoint regions, or boundaries of the 
AR, are characterized by a blueshift of ^20 km s~^ and an enhancement of the line width. 
The deviation from a single Gaussian profile, which is quantified by the enhancement of Xr^ 
is also clearly seen in the loop footpoint regions. The Si X 258.37A and Fe xiii 202.04A 
lines reveal basically the same structures in the maps of Gaussian parameters. The cool line 
Si VII 275. 35A shows prominent red shifts in the fan regions (not shown here), a phenomenon 
found also by I Warren et al.l ( 120111 ) . Such redshifts could be the downfiows a fter cooling of 



the high-speed upfiows and evidence for the mass circulation in coronal loops ( iMarsch et al. 
I2OO8,). 



4.1. RB asymmetry analysis 

Similarly to before, we applied the RBp technique to each spectrum and calculated the 
average relative intensity in the velocity interval of 60-110 km s~^ from the RB asymmetry 
profile. Figure [7](H) shows the spatial distribution of this average RB asymmetry. It is 
clear that the blueward asymmetry is most prominent at the two edges of the AR, generally 
coincident with the enhancement of blueshift, line width, and Xr- However, the blueward 
asymmetry is clearly not only present at the AR edges, but also existing in other plage/loop 
footpoint regions. For example, the tongue-like patch (x=— 230'',y=— 430^') of blueward 
asymmetries in Figure [7](H) seems to be coincident with moss and it is not a typical outfiow 
region. The fact is that simple magnetic field structures which are aligned with the line of 
sight are often not present in AR cores so that the secondary component there is often not 
revealed as clear blueward asymmetries of line profiles. An obvious correlation between the 
Doppler shift/Gaussian width derived from single Gaussian fit and peak relative intensity 
derived from RBp asymmetry profiles is found in Figure [HI suggesting that the blueshift 
and enhanced line width are caused by the high-speed secondary emission component and 
that the single Gaussian fit can not accurately refiect the real physical process here. In the 
following, only those pixels with obvious blueward asymmetry (the average RB asymmetry 
in 60-110 km s~^ smaller than -0.05) and significant signal to noise ratio (larger than 5) were 
selected. Clearly, such a cutoff excludes the analysis of spectra in the AR core and edges 
where the magnetic field lines make large angles relative to the line of sight. The peak relative 
intensity, velocity, and 1/e width derived from RBp asymmetry profiles of Fe XIV 274.20A 
at these pixels are shown in Figure [71(E)- (G). 
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The distributions of these three parameters are shown as the red histograms in Fig- 
ure [lT]( A)- (C) and they are very similar to those in Figure [6l The only obvious difference 
is perhaps the slightly larger peak velocity and Gaussian width. In the 2007 January 18 
observation, the velocity and width peak at ^75 km s~^ and ^50 km s~^, respectively. But 
in the 2010 September 16 observation, the velocity peaks at ^85 km s~^ and the width peaks 
at ^60 km s~^. The major part of the velocity distribution is still on the left side of 100 
km s~^ in the 2010 September 16 observation, which is no t consistent with the scenario of 
a considerably broader secondary component ( |Peterll2010l ). The scenario of a considerably 
narrow secondary component can also be excluded since the RBp width is almost always 
larger than 35 km s~^. Similar to the 2007 January 18 observation, the RBp velocity in this 
recent observation also seldom reaches higher than 150 km s~^, indicating that the line of 
sight velocity of the secondary component is unlikely to be as high as 200 km s~^ in this 
observation. So with the RB asymmetry analysis (particularly RBp), we can conclude that 
the secondary emission component, which is usually less than 20% and often a few percent 
of the primary one in intensity, has a speed around ^80 km s~^ and a width comparable to 
the primary component. 



4.2. Double Gaussian fit 

The same algorithm of RBp guided double Gaussian fit was adopted to the asymmetric 
profiles obtained in the 2010 September 16 observation. Several examples of the observed 
and fitted line profiles are presented in Figure M In Figure M^A) we show profiles of four 
emission lines averaged over region 1 which is marked in Figure [71 Blueward asymmetries 
and deviations from single Gaussian profiles are clearly present in all of the four lines. The 
double Gaussian algorithm yields a reasonably good fit to the observed spectra. This spatial 
average is necessary because individual profiles of the Si Vll 275. 35A and Si x 258. 37A 
lines are too noisy to allow reliable Gaussian fits and RB asymmetry analysis in the AR 
boundaries. 

We find that the relative intensity and velocity of the secondary component stay rela- 



of multi-thermal hi^h-speed upfiows, as suRRested by our 


previous 


work ( 


De Pontieu et al. 


2009: De Pontieu & Mcintosh 2010: 


De Pontieu et al. 2011: 


Mcintosh & De Pontieull2009a b: 


Mcintosh et al. 


2010a, b: 


Tian et al. 


2011a 


b: Martinez- Svkor a et al. 


boila 


). The clear pres- 



ence of blueward asymn i etry of the cool Si Vll 275. 35A line is consistent with the finding of 
Mcintosh fc De Pontieul ( j2009al ). Compared to other lines, the Si Vll 275. 35A line reveals 
a smaller velocity (^65 km s~^) of the secondary component. One possible reason could 
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be that the Si Vll 275. 35A hne reflects emission from much cooler plasmas so that its line 
profiles are complicated by cooling downflows (see also Ugarte-Urra fc WarrenI 2011 ). 



Fe XIV 274. 20A line profiles at four single pixels (without spatial average) are presented 
in Figure [9](B)-(E). The secondary component resolved by the double Gaussian fit has dif- 
ferent relative intensity and velocity. The magnitude of velocity is likely to be related to the 
viewing angle of the upfiows. If the magnetic field lines associated with the upfiows make 
a smaller angle with respect to the line of sight, the observed velocity should be closer to 
the real velocity. The location of the profile presented in Figure [9](B) is in a weak emission 
region surrounded by hot loops. This weak emission region might be related to less compact 
(larger) loops or open field lines which are usually more vertical at lower heights compared 
to compact and smaller loops. This may explain the large value of velocity (^120 km s~^). 
The Gaussian widths of the two components are generally comparable in all these four cases. 

The KBd profiles were then calculated by taking the velocity of the primary component 
as the line centroid. In Figure [9] we also show the three RB asymmetry profiles (RB^-, RBp, 
RBd) for the corresponding spectra. The parameters of the secondary component obtained 
through RBp and RBp) techniques are usually close to those obtained by using the double 
Gaussian fit. 

The RB p guided double Gaussian fit and RB ^ analysis were then applied to all spectra 
of Fe XIV 274. 20A with obvious blueward asymmetries. Figure [TOl shows the spatial distri- 
butions of the three Gaussian parameters of both components as well as the intensity ratio 
and velocity difference of the two components. Again, we found that the primary compo- 
nent is usually slightly blue shifted by ^10 km s~^ in this observation. The velocity of the 
secondary component is blue shifted by much larger values. The maps of parameters derived 
from RBp) asymmetry profiles are very similar to those derived from RBp, which have been 
presented in Figure [3 and thus are not shown here. Structures on the maps of RBp peak 
relative intensity, velocity, and width shown in Figure [7] generally coincide with those on the 
maps of intensity ratio, velocity difference, and secondary component width as presented in 
Figure [TOl indicating that both the RBp technique and RBp guided double Gaussian fit are 
able to extract the parameters of the secondary component. 

Figure [TT](A)-(C) shows the distributions of the intensity ratio (relative intensity), ve- 
locity difference (velocity), and Gaussian width of the secondary component as derived from 
double Gaussian fit and RB asymmetry analysis (RBp,RBp)). For each of the parameters, 
the three histograms are very similar. All these histograms show no big difference from those 
presented in Figure [6l except the slightly larger peak velocity and width. The velocity of the 
secondary component is usually in the ra nge of 50-12 km s~^ in this observation, which is 



about twice as high as that suggested by IPeterl (120101 ) . 
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The distribution of the Xr ratio between the double and single Gaussian fit in Fig- 
ure [TTlfD) is also very similar to that in Figure [6](D). The ratio is usually less than 0.5, 
indicating a significant improvement of the fitting by using the double Gaussian fit instead 
of the single Gaussian fit. 



4.3. AIA observations 

During the EIS scan, the AIA instrument took full-disk images of the Sun in several 
EUV channels at a cadence of 12 s. We selected images in the passbands of 171A, 193A and 
21lAfor our analysis. The dominant emission of these passbands has a temperature around 
0.8 MK, 1.3 MK, and 2 MK, respectively. There were a few missing frames in each passband 
so that an linear interpolation was applied to the data. The pixel size is about 0.&\ 

We extracted a sub-region enclosing the AR observed by EIS and co-aligned different 
imase frames. A movie of AIA 193A images and the running diflFerence images is available 
online (AIA193.mpeg). Figure [T2l( A) shows an image of AIA 193A at 12:06. The rectangle 
indicates the region scanned by EIS. The cross correlation technique was used to coalign 
AIA 193A and EIS Fe xiv 274.20A. 

Upward propagating disturbances (PDs) are clearly present in the AIA movie. These 
disturbances are mainly found at the two boundaries of the AR and coincide with obvious 
blueward asymmetries as seen from Figure [71(H), suggesting a connection of the two. The 
speed of these disturbances can be estimated by placing a virtual slit along the propagation 
direction of the disturbances and calculating the slope of a bright strip in the space-time 



(S-T) plot fe.g. JSakao et al.ll2007l : Mcintosh et al.ll2010al : iTian et al.ll2011bl ). Figure[l2](B) 



(G) show the S-T plots of the original and detrended intensities for the dotted line shown in 
Figure [T2l( A). The detrended intensities were used to better reveal the faint outfiow signatures 
and they were obtained by first subtracting a 8-minute running average from the original 
intensity time series and then normalized to the running average at each location of the slit. 
One example of PD is indicated by the inclined dashed line. The speed of this PD, the slope 
of the dashed line, is estimated to be 91 km s~^. 

As the PDs are usually present at different passbands of AIA, in principle we should 
be able to compare the speed of a certain PD at different temperatures using AIA data. 
For a certain virtual slit, we often find similar patterns of the S-T plots (both the orig- 
inal and detrended intensities) in the 193 A passband and the other two passbands (par- 
ticularly the 211 A passband). In most cases we did not find an obvious increasing trend 
in the speed with increasing temperature, which seems to support the interpretation of 
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PDs as multi-thermal hiRh-speed upflows (e.R., 



2009a: 



De Moortel et al 



2000, 



De Pontieu fc Mclntoshll2010l : iTian et al 
2OO2I: 



Robbrecht et al.ll2001 



Sakao et all 120071 : iMcIntosh fc De Pontieu 



201 laf) rather than slow mode waves (e.g., 



^mg et al.ll2003l:lMcEwan fc De Moortel 



Verwichte et al.T201ol : Istenborg et al.l 



20061 : 



Marsh et al. 



I2OO9I . 



2011 



Wang et al. 



2009a 



20111). However, a recent investigation suggests that non-dominant cool ions s uch as V may 

contribute significantly to the total emission of the 193A and 211 A passbands (IMartmez-Sykora et al. 
2011bl ). But cool ions such as O V would not be expected to show obvious emission over a 
large height range. For example, PDs are also seen off-limb in a coronal hole up to heights 
of 50-100 Mm with AIA ( iTian et al.ll2011bl ). i.e., heights where emission from cool lines does 
not typically occur at sufficient strength. This renders the impact of non-dominant cool 
ions less likely. However, we can not rule out the possibility of slow waves. Perhaps both 
slow waves and upfiows are existing in our observations. In fact waves can be exited by 
high-speed flows. But our EIS result that lines with different formation temperatures basi- 
cally show similar velocities seems to suggest that upflows dominate the emission. We note 
that for some virtual slits the S-T plots show some differences in different AIA passbands, 
especially between the 171 A passband and the ot her two passbands. T he emission in the 
171 A passband is dominated by Fe ix 171.107A ( lO^Dwyer et al.ll2010l ). which is formed 
in the upper transition region and thus may have a smaller scale height compared to the 
hotter emission in the other two passbands. Moreover, due to the different sensitivity, signal 
to noise ratio, and emission contrast in different passbands we may not always be able to 
clearly identify every outflow event in all the three passbands from the S-T plots. The 211 A 
passband samples plasma with a temperature much closer to the formation temperature of 
Fe XIV 274.20A. However, the 211 A images have a lower signal to noise ratio and are thus 
noisier compared to the 193 A images. Since the S-T plots in the 211 A passband are highly 
similar to those in the 193 A passband, in the following we concentrate on the 193 A images. 

We identifled 39 well-isolated disturbances propagating along the fan-like structure at 
the left boundary of the AR and calculated their propagating speeds in the 193 A passband. 
A speed distribution was then obtained. A comparison of this speed distribution with the 
distribution of the velocities derived from EIS Fe xill 202. 04 A line proflles (through double 
Gaussian flt, RBp, and RB^ analysis) is presented in Figure [T3l( A) & (B). Here only line pro- 
flles with obvious blueward asymmetries in the fan-like structure at the left boundary of the 



AR were used. By using the P otential Field Source Surface (PFSS) model ( ISchatten et al 



19691 : ISchrijver fc DeRosall2003l ). the inclination angles of the magnetic fleld lines in the fan- 
like structure at the left boundary of the AR were estimated to be around 55° with respect 
to the line of sight. Disturbances propagating at a speed of V along these magnetic fleld 
lines should exhibit a plane-of-sky speed of Vsin55° and a line-of-sight speed of Vcos55''. 
The ratio of the two speed components is 1.43. Since the speeds derived from S-T plots 
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of AIA images are the real speeds projected onto the plane of sky and the speeds derived 
from EIS line profiles are the line-of-sight component of the real speeds, the ratio of the two 
should be around 1.43. From Figure [T3l we find that this is indeed the case. The speeds 
derived by using AIA data are systematically larger than those derived by using EIS data. 
The ratio of the two average speeds is around 135/85=1.58, which is close to 1.43. For com- 
parison, the RBp velocity histogram is multiplied by 1.58 and overplotted as the dotted line 
in Figure [T3l(B). The small difference between 1.58 and 1.43 could easily be caused by the 
fact that the PFSS model is not expected to give a perfect representation of w hat viewing 



ang:les the real Sun presents. W e think that in the future 3-D reconstruction (e.g. JFeng et al 
20071 : lAschwanden et al.ll2008l ) could be considered as an alternate method to help relate the 
velocities derived from imaging and spectroscopic observations. In Figure [13] we can also 
compare histograms of the relative intensities derived from EIS Fe xill 202. 04 A line profiles 
(in panel (C), through double Gaussian fit, RBp, and RBd analysis) and those derived from 
S-T plots (D). The latter was obtained through the following process: we first smoothed 
the original time series with a 1-hour kernel. After that we subtracted the smoothed time 
series from the original one and obtained a time series of intensity change. The relative 
intensity was then obtained by normalising the maximum intensity change to the average 
intensity at each slit location. From Figure [T3l( D) we can see that the relative intensity is 
usually a few percent but can also be larger than 10%, which is consistent with magnitude of 
the relative intensity of the secondary component revealed by RB asymmetry analysis and 
double Gaussian fit. For a better comparison, the RBp relative intensity histogram is over- 
plotted as the dotted line in Figure [T3l(D). Note that the lack of small values of the relative 
intensity derived from EIS line profiles are due to the fact that we only selected profiles with 
significant blueward asymmetries. These results provide further support to the argument 
that the blueward asymmetries o f coronal line profiles are caused by the upward pr opagating 



disturbances at AR boundaries ( IMcIntosh fc De Pontieull2009al : iTian et al.ll2011al ). 



As mentioned bv iDe Moortell ( 120091 ). although there are striking similarities, it is not 
clear how the outfiows in ferred from EIS lin e profiles and the PDs seen in EUV imaging 



observations are related. iDe Moortell ( 120091 ) pointed out that the blue shifts (from single 



Gaussian fit) of EIS lines are only of the order of 20-50 km s~^ and no periodicity has been 
reported. While the PDs in imaging observations have a speed of ^100 km s~^ and they are 
often quasi-periodic. Now we understand that the EIS line profile consists of two components 
and that the secondary component has a velocity of ^100 km s~^. And we have already 
demonstrated that the secondary component is a lso quasi-periodically enhanced or weakened 
( IDe Pontieu fc Mclntoshll2010l : iTian et al.ll2011al ) . The consistency between properties of the 
PDs observed by AIA and those of the secondary component revealed by RB asymmetry 
analysis and double Gaussian fit of EIS line profiles suggests a close connection of the two. 
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It is thus natural to suggest that the PDs are responsible for the blueward asymmetries of 
line profiles. The non-thermal width of the secondary component is probably caused by both 
the superpos ition of multiple unresolved upflows and the Alfvenic motions associated with 
the upfiows dMcIntosh fc De PontieulboOQal : IPoUa fc Zhukovlboilh . 



We have to mention that in Figure [T3l we only present the results for the fan-like structure 
at the northeast boundary of the AR because of the simple magnetic field structure there. 
From the online movie (AIA193.mpeg) we can see that outward propagating disturbances 
are also clearly present in the fan-like structure to the southwest of the AR. We produced 
several S-T plots and found generally similar S-T patterns in diflFerent passbands. However, 
the footpoints of these fan loops (around x=— 110^', y=— 45C) are mixed together and the 
higher parts of diflFerent loops are not well-separated, making it difficult to identify a well- 
isolated propagating disturbance in the S-T plots. In other parts of the AR, magnetic field 
structures are either too complex or only slightly inclined with respect to the line of sight so 
that reliable speeds may not be obtained through S-T plots. 



5. Conclusion 

In conclusion, the presence of asymmetric coronal line profiles in loop footpoint regions 
suggests that at least two emission components are present: an almost stationary primary 
component and a high-speed secondary component. The secondary emission component 
might play an important role in supplying the corona with hot plasma. 

We have generated artificial spectra composed of two Gaussian components to test the 
ability of extracting the Gaussian parameters (intensity, velocity, and width) of the sec- 
ondary component through the originally defined (use the single Gaussian fit to determine 
the line centroid) and a slightly modified (use the spectral position corresponding to the 
peak intensity as the line centroid) technique of the RB asymmetry analysis, and find that 
this technique, especially the modified one, can provide a relatively accurate estimate of the 
parameters of the secondary component if the velocity is larger and the width is not consid- 
erably larger than the primary component width. We have applied the modified technique 
to the spectra obtained by the EUV Imaging Spectrometer (EIS) onboard Hinode in two 
active regions and found that the distributions of the velocity /width determined from the 
RB asymmetry analysis peak around 80 km s~^/60 km s~^. A comparison of such distribu- 
tions with the values of velocity/ width calculated from the artificial spectra suggests that 
the secondary component cannot be very broad or narrow and that the widths of the two 
components are comparable. The velocity of the secondary component is usually within the 
range of 50-150 km s~^ and seldom reaches as high as 200 km s~^ or as low as 40 km s~^ in 
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our EIS observations. The relative intensity of the secondary component with respect to the 
primary one is often a few percent but can sometimes reach 30%. Our conclusions are very 



different from those of iPeterl (120101 ). who claimed that the secondary component contributes 
some 10% to 20% to the total emission, is about twice as broad as the primary component 
and blueshifted by up to 50 km s~^. The velocity range of the secondary com ponent we de- 
rived here is roughly consistent with, yet slightly smaller than that derived by Br vans et al 



(2010). However, we find that the very high speed of 200 km s ^ mentioned by 



Bryans et al 



( 120101 ) seems to be not very likely. 



Using the parameters determined from the RB asymmetry analysis and single Gaussian 
fit as initial values for several free parameters, we have performed a double Gaussian fit to 
the spectra that were observed to have obvious blueward asymmetries. The RB asymmetry 
analysis has been performed for the third time, but then using the centroid of the sec- 
ondary component determined from the double Gaussian fit as the line centroid. The double 
Gaussian fit and the RB asymmetry analyses have given basically consistent results. The 
properties of the secondary component stay relatively stable across the temperature range 
of log(T/K)=5.8-6.3. The double Gaussian fit also shows that the velocity of the primary 
component is often blue shifted by ^10 km s~^. 

We have also used imaging observation simultaneously performed by the AIA instrument 
onboard SDO and demonstrated that the propagating disturbances coincide with obvious 
blueward asymmetries of line profiles, show an intensity change similar to the relative in- 
tensity of the secondary component of line profiles and an average velocity consistent with 
that derived from emission line profiles. Such results suggest that the upward propagating 
disturbances are plasma upfiows responsible for the blueward asymmetries of line profiles. 



SDO is the first mission of NASA's Living With a Star (LWS) Program. EIS is an 
instrument onboard Hinode^ a Japanese mission developed and launched by ISAS/JAXA, 
with NAOJ as domestic partner and NASA and STFC (UK) as international partners. It is 
operated by these agencies in cooperation with ESA and NSC (Norway). Scott W. Mcintosh 
is supported by NASA (NNX08AL22G, NNX08BA99G) and NSF (ATM-0541567, ATM- 
0925177). Bart De Pontieu is supported by NASA grants NNX08AL22G and NNX08BA99G. 
Hui Tian is supported by the ASP Postdoctoral Fellowship Program of NCAR. The National 
Center for Atmospheric Research is sponsored by the National Science Foundation. Hui Tian 
thanks Marc DeRosa for his assistance in using the PFSS package. We appreciate the efforts 
of the referee to improve the manuscript. 
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A. A comparison between RBp and RB^^ results 



As we mentioned in Section 2, both the RBp and RB^- techniques can provide a relatively 
accurate estimate of the Gaussian parameters of the secondary component when the offset 
velocity is larger than the Gaussian width of the primary component. We have also found that 
the modified RB technique, RBp, has a better ability to accurately quantify the properties 
of the secondary component as compared to the originally defined RB^- technique. Figure [T4l 
shows a comparison between the RBp and RB^- results for the 2010 September 16 observation. 
We can see that the spatial structures in all of the four parameters are very similar. The 
most prominent diflFerence is the smaller values of relative intensity and width, and larger 
upfiow velocity for RB^'. Such results are fully consistent with our forward modeling results 
in Figure [1] and imply that in our previous work the relative intensity and width might be 
underestimated while the upfiow velocity was likely to be slightly overestimated. 



B. A comparison between EIS Fe xiv 274. 20A and Fe xiv 264. 78A line 

moments and profile asymmetries 

We have mentioned that the strong Fe xiv 274. 20A line is blended with the weak 
Si VII 274. 18 A lin e and that the latter can safely be ignored in active region conditions 



( lYoung et al.l 120071 ). The very close wavelengths (less than 1 spectral pixel) of the two also 
suggest that the blend is very unlikely to have an important impact on the high-speed (of 
the order of 100 km s~^) upflow. The Fe xiv 264. 78A line spectra were not obtained in the 
2007 January 18 observation so that we can not provide a direct comparison of it with the 
Fe XIV 274. 20A line. In 2010 September 16 both Fe xiv lines were used in the observation 
and Figure [T5l shows the spatial distributions of the line moments and profile asymmetries. 
Note that the spectra (especially for the Fe xiv 264. 78 A line) in the upper left corner, 
the higher part of the fan, have low signal to noise ratio and thus the line moments and 
profile asymmetries are noisy there. The very similar behavior of line moments and profile 
asymmetries between Fe xiv 274. 20A and the weaker but clean Fe xiv 264. 78A line provides 
further confidence that the blend of Si Vll 274. 18A line can be safely ignored in our analysis. 



C. EIS observation on 2007 December 12 



We also applied the techniques of RB asymmetry analysis and double Gaussian fit 
to an EIS observat ion of AR NOAA 10978 from 11:43 to 17:03 on 2007 December 12. 
Bryans et al.l (120101 ) found that the secondary emission component in this AR has veloc- 
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ities around 1 10 km s ^ and, as they stated in their paper, often as high as 200 km s 



Bryans et al.l (120101 ) used the Fe xil 195. 12 A hne and they mentioned that the blend of this 



l ine, Fe xil 195. 18 A, should not significantly a ffect the Fe xil 195. 12A line in the low-density 



( IDoschek et al.ll2008l : iBrooks fc Warrenll201ll ) outflow region (AR edges). However, to avoid 



the complexity introduced by this blend we decided to turn to other clean and strong Fe xil 
lines. We found that in the AR edges some other lines at the far blue wing of the Fe xil 
192.39A line can be relatively strong and thus aflFect the results of RB asymmetry analysis 
and double Gaussian flt. So here we mainly present the results of the Fe xil 193. 5lA line, 
which is clean and about 2/3 as strong as Fe xil 195. 12A. Since the line is very strong, there 
is no need to perform a spatial average of the line proflles. 

Figure [H] shows the spatial distributions of the single Gaussian parameters and proflle 
asymmetry. Similar to AR 10938 and AR 11106, the two loop foot point regions are char- 
acterized by an enhancement in the blue shift, line width, and blue ward proflle asymmetry, 
and there is an obvious correlation between the Doppler shift/Gaussian width and proflle 
asymmetry. Note that the white strip in each panel mark locations where the observed data 
is lost. 

Two examples of line proflles and the corresponding RB asymmetry proflles of the 
Fe XII 193.51A and Si vii 275.35A lines are shown in Figure \17\ The proflles were averaged 
over 3 pixels center ed at the locations marked in Figure [T6l(A) in both spatial dimensions. 



Bryans et al.l (2010) mentioned that cool lines like Si Vll 275. 35A do not have any signiflcant 
asymmetries and that their proflles can be accurately represented by a single Gaussian. We 
found that this is not always the case. Similar to what we found in AR 11106, at some 
locations the Si Vll 275. 35A line proflles clearly exhibit obvious asymmetries. As shown in 
Figure [TTl usually the primary component is redshifted by ^15 km s~^ and the velocity 
difference between the two components is smaller than that of Fe xil 193. 51A, which might 
be due to the complexity introduced by cooling downflows. 

Similar to what we did for AR 10938 and AR 11106, the techniques of RB asymmetry 
analysis and RBp guided double Gaussian flt were applied to each spectral proflle with an 
obvious blueward asymmetry (the average RB asymmetry in 60-110 km s~^ smaller than 
-0.05). Spatial distributions of the double Gaussian parameters are presented in Figure [T8l 
Spatial distributions of the Gaussian parameters derived from the RBp and RB^ analyses 
are similar to those of the secondary component shown in Figure [IHl and thus are not shown 
here. 

The distributions of the intensity ratio (relative intensity), velocity difference (velocity), 
and Gaussian width of the secondary component as derived from double Gaussian flt and 
RB asymmetry analysis (RBp,RBp)), as shown in Figure [T9l( A)- (C). reveal no big difference 
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from those in Figure [TTl However, the velocity distributions seem to peak at ~95 km s"-*^, a 
value larger tha n those in AR 10938 (^75 km s~-^) and AR 11106 (^85 km s~-^) and close to 
that derived by Bryans et al. ( 201ol ) (^110 km s~M. Th e velocity seldom reaches as high as 
200 km s~^. In fact from Figure 7 of iBryans et al.l (120101 ) one can conclude that the velocity 
of the secondary component is usually in the range of 60-170 km s~^ and that only a very 
minor portion of the velocities exceeds 200 km s~^. 

As seen from Figure [T9l(D). the Xr ratio between the double and single Gaussian fit 
is generally smaller than 1, indicating an improvement of the fitting by using the double 
Gaussian fit instead of the single Gaussian fit. 
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Fig. 3. — Scatter plots of the relationship between Doppler shift/Gaussian width derived 
from single Gaussian fit and peak relative intensity derived from RB asymmetry analysis 
(RBp) for Fe xiv 274.20A. 
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(A) FeXIV274.20 line profile FeXV284.1 6 line profile (B) FeXIV274.20 line profile FeXV284.1 6 line profile 
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Fig. 4. — (A) RB asymmetry profiles of the Fe xiv 274. 20A and Fe XV 284. 16A line profiles 
averaged in region 1 marked in Figure [2]( A) . Top: the observed spectra and measurement 
errors are shown as the diamonds and error bars, respectively. The dotted lines are single 
Gaussian fits. The two dashed lines in each panel represent the two Gaussian components 
and the solid line is the sum of the two components. The velocity (v) and Gaussian width 
(w) derived from the single (SGF) and double (^^2^^ for the two components) Gaussian 
fits are shown in each panel. Also shown is the intensity ratio of the secondary component 
to the primary one (i2/ii)- Bottom: the solid, dashed, and dotted lines represent RB profiles 
for RB^-, RBp, and RBp), respectively. The peak relative intensity (i), velocity (v), and 1/e 
width (w) are shown in each panel. (B) Same as (A) but for the line profiles averaged in 
region 2 marked in Figure [2]( A) . 
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(A) r' Peak Int. (arbitrary unit) (^) r' Doppler Shift (km/s) (^) 1 ^' Gaussian Width (km/s) (P) Intensity Ratio 
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(E) Peak Int. (arbitrary unit) (F) 2"' Doppler Shift (km/s) (G) 2"' Gaussian Width (km/s) (H) Velocity Difference (km/s) 
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Fig. 5. — Parameters of the primary (A-C) and secondary (E-G) components as derived 
from the double Gaussian fit for Fe xiv 274. 20A. The intensity ratio and velocity difference 
of the two components are shown in (D) and (H), respectively. 
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Fig. 6. — Histograms of the intensity ratio (A), velocity difference (B), and Gaussian 
width of the secondary component (C), as derived from double Gaussian fit (black) and 
RB asymmetry analysis (red/blue for RBp/RB^)) for Fe xiv 274. 20A. Panel (D) shows the 
histograms of the Xr ratio between the double and single Gaussian fit. 
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(A) Peak Intensity (arbitrary unit) (B) Doppler Sliift (km/s) (C) Gaussian Widtli (km/s) (D) of Single Gaussian Fit 
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Fig. 7. — Same as Figure [2] but for the observation on 2010 September 16. The square in 
(A) mark the location where profiles of several emission lines are averaged and presented in 
Figure [9]( A) . The four pluses mark locations of the Fe xiv 274. 20A line profiles presented in 
Figure [9KB)- (E) 




Fig. 8. — Same as Figure [3] but for the observation on 2010 September 16. 
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Fig. 9. — (A) RB asymmetry profiles (bottom) of four emission line profiles (top) aver- 
aged in the square marked in Figure [7]( A) . (B)-(E) RB asymmetry profiles (bottom) of the 
Fe XIV 274. 20A line profiles (top) at locations 2-5 marked in Figure [71(A) . The line styles 
and denotations of parameters are the same as in Figure [H 
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Fig. 10. — Same as Figure [5] but for the observation on 2010 September 16. 
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Fig. 11. — Same as Figure [6] but for the observation on 2010 September 16. 
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Fig. 12.^ (A) The image of AIA 193A at 12:06 on 2010 September 16. The rectangle 
outhnes the region scanned by EIS. (B) A space-time plot of the AIA 193 A intensity for the 
dotted line shown in (A). (C) Same as (B) but showing the detrended intensity. (D) & (E): 
Same as (B) & (C) but for AIA 17lA. (F) & (G): Same as (B) & (C) but for AIA 21lA. The 
inclined dashed line indicates one propagating disturbance with a traveling speed of about 
91 km s"-*^. A movie associated with (A) is available on line (AIA193.mpeg). 
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Fig. 13. — Histograms of the velocity and relative intensity derived from the EIS 
Fe XIV 274.20A line profiles (A, C) and AIA 193A S-T plots (B, D) in the fan-like region 
at the left boundary of NOAA AR 11106. In (A) & (C) the black, red and blue histograms 
represent results of the double Gaussian fit, RBp, and RBd analysis, respectively. For com- 
parison, the RBp velocity histogram is multiplied by 1.58 and overplotted as the dotted line 
in (B), and the RBp relative intensity histogram is overplotted in (D) as the dotted line. 
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Fig. 14. — Spatial distributions of the parameters derived from RBp (A-D) and RB5' (E-H) 
asymmetry analysis for Fe XIV 274. 20A in the observation on 2010 September 16. 
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Fig. 15. — Spatial distributions of the line moments (single Gaussian parameters) and profile 
asymmetries (average relative intensity in the velocity interval of 60-110 km s~^, as obtained 
from the RBp profiles) for Fe xiv 274. 20A (A-D) and Fe xiv 264. 78A (E-H) in the observation 
on 2010 September 16. 
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Fig. 16. — Spatial distributions of single Gaussian parameters and profile asymmetries (av- 
erage relative intensity in the velocity interval of 60-110 km s~^, as obtained from the RBp 
profiles) for Fe xil 193. 51A in the observation on 2007 December 12. The two pluses mark 
locations of the line profiles presented in Figure [T71( A) &(B), respectively. 
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Fig. 17. — (A) RB asymmetry profiles of the Fe xil 193. 51A and Si Vll 275. 35A line profiles 
at location 1 marked in Figure [T6l(A). The line styles and denotations of parameters are 
the same as in Figure [H (B) Same as (A) but for the line profiles at location 2 marked in 
Figure [ISK A). 
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Fig. 18. — Same as Figure [5] but for the observation on 2007 December 12. 




Fig. 19. — Same as Figure [6] but for the observation on 2007 December 12. 



